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Abstract 

We present a simple technique that allows to generate Feynman 
diagrams for vector models with interactions of order 2n and simi- 
lar models (Gross-Neveu, Thirring model) using a bootstrap equation 
that uses only the free field value of the energy as an input. The 
method allows to find the diagrams to, in principle, arbitrarily high 
order and applies to both energy and correlation functions. It auto- 
matically generates the correct symmetry factor (as a function of the 
number of components of the field) and the correct sign for any dia- 
gram in the case of fermion loops. We briefly discuss the possibility 
of treating QED as a Thirring model with non-local interaction. 
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1 Introduction 



In field theory the calculation of physical quantities is commonly given as 
a sum of Feynman diagrams. Each diagram is calculated with a certain 
weight, dependent on the number of symmetries of the diagram. A rule to 
find this symmetry factor is given in 0, [|. Furthermore expressions were 
found for the total number of diagrams that contribute in any order to various 
quantities of interest 0, H . In order to find the number and type of diagrams 
that contribute to a physical quantity, numerous programs were developed 
|@, 0]. In this article we want to present a technique for generating Feynman 
diagrams that is a generalisation of a similar technique previously used in 
the study of the vector model with quartic interaction in 0+0 dimensions [p]]. 
The results obtained in this paper are valid for any number of dimensions. 
We will find all diagrams that contribute to the energy and the correlation 
functions (to, in principle, arbitrarily high order) with the correct symmetry 
factor. This factor can be a function of the number of components of the 
field. The technique is essentially different from the ones that are usually 
applied tot generate diagrams, since it does not involve any combinatorics. 
It is a recursive method that requires only the free field value of the energy 
as an input. 

In m it was observed that a particular relation between the partial deriva- 
tives of the energy of the vectormodel could be used to obtain a perturbation 
series for the energy. In the simple case of the 0+0 dimensional vector model 
it also allowed to find explicit expressions for the coefficients of the large N 
expansion of the model. A general form for this coefficients could be found 
and proven to all orders. Unfortunately, a straightforward application of the 
method used in 0+0 dimensions to higher dimensions is not possible. How- 
ever, if we slightly adapt the original model a similar technique can be used 
in higher dimensions. We will consider here the vector model with quar- 
tic self-interaction. Instead of considering a model with constant mass, we 
will allow the mass to depend on the coordinates (the mass is therefore an 
external field). The model is defined by the partition function: 

Z = e- E = Jv ( f ) exp-^Jdxd^d II (P + m(x) ( f ) 2 (x)-gJ( ( j ) 2 ) 2 , (1) 
where </> is a N-component field. If we now consider the first and second 
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functional derivative of E with respect to m(x) we find: 

9E 1 Alt \ 

- < (f) (x) > 



dm(x) 2 

and: 

S ' E 1 - ,2, U2 - , , 9E N2 



where < ... > denotes the expectation value of an operator. Integrating the 
last equation with respect to x we find: 



4 / (te^r) 2 - T^frr) (4) 



We rewrite this equation in the form that will be of interest to us: 

dE A r., dE x2 d 2 E 

= 4 / (( ) 

dg J dm(x) dm(x) 

This is the generalisation of the formula found in [l]. It is clear how this 
formula can be used to generate a perturbation series for the vectormodel. If 
we start with the zeroth order approximation of the energy (i.e.: the value of 
the energy when there is no interaction), formula (|]) gives us the derivative of 
the energy (with respect to the coupling constant) to the same order. Upon 
integration, we therefore find the energy to first order. This again can be used 
in (f|) to find the energy to second order. By continuing this way we generate 
a perturbation series solely from the knowledge of the free theory and the 
fact that the energy obeys @). Of course, since the perturbation series is 
known to be given by a sum of Feynman-diagrams it is not unexpected to 
find that we will simply generate the diagrams to any desired order, with the 
exact symmetry factors. 



2 Generating Feynman diagrams 

In order to be able to carry out the procedure described in the previous 
section, we must know the value of the energy when there is no interaction. 
Without any difficulty this is found to be: 

E(g = 0) = — Tr log(-d 2 + m) (5) 
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(remember that m is a field, rather than a constant) If we denote by {\x >} a 
complete set of eigenstates of the position operator, then (|5|) can be written 
as: 

N f m 
E{g = U) = - J dy<y\\og{l-—)\y> (6) 

For later convenience we have substracted a quantity that is independent of 
m. This will not change the result. We will also denote — d~ 2 as G. The 
expression (f|) is defined by its series expansion, and so is the functional 
derivative with respect to m(x). If we consider a general term in the series 
expansion, the derivative will work successively on each of the m's that occur 
in the term. We find that the derivative of such a general term is given by: 

3 n ~ 1 

<y\{mG) n \y>=Y,<y\{™Gy\x><x\G{mG) n ~ 1 - 1 \y> (7) 

j=0 



dm(x 

And when we integrate with respect to y: 
d 

dm(x) 
We therefore find: 



J dy < y\(mG) n \y >= n< x\G(mG) n - l \x > (8) 



-— - < x\G- —\x > (9) 



dm(x) 2 1 + mG 

We will also need higher derivatives of the energy. Let us therefore con- 
sider the following derivative: 

dm(x) < + > dm(x) < ^^J^/ mG) \z > (10) 

Again using (|?|), we find that this is equal to: 

oo n 

~ E E < y\G(-mG) j \x >< x\G(-mG) n - J \z > (11) 

n=0jr'=0 

By rearranging the double sum we can write this in the form: 

- < v\G ^— -\x >< x\G —=\z > (12) 
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This formula has a simple interpretation. Since < y\G 1+mG \z > is nothing 
but the propagator of the field between points y and z, taking the derivative 
with respect to m(x) is seen to double this propagator by inserting an extra 
point x between the two endpoints (apart from a change of sign). We can 
denote this graphicly as: 



We will now illustrate the techique by generating the first orders in the 
series expansion of the energy. Starting from the lowest order approximation 
(i.e. the value of the energy without interaction), we find from the formulas 
given above: 

dE N 1 

— — - = — <x\G- -\x>+0(g) (13) 

dm(x) 2 1 1 + mG 1 KJ v ' 

d2E N ^ 1 , 1 , 

— < x\G — \x >< x\G — \x>+0(g) (14) 



dm(x) 2 2 1 + mG 1 + mG 



Using (p73|) and in (01) and integrating with respect to g (taking into 



account the correct integration constant), we find: 

E = —Tr\og(l + mG) + N(N + 2)g [ dx < x\G l -—:\x > 2 +0(g 2 ) (15) 

2 J 1 + mG 

As long as m is a function of x the integral can not be carried out. Unfortu- 
nately, we can not put m equal to a constant at this point of the calculation, 
because we need the full expression to find the higher order terms. It is clear 
however that the first order expression we have found corresponds to the 
diagram: 



E=f()+N(N + 2)g{_XJ +0(g 2 ) 



as can be seen by drawing the propagators that occur in the term or 
by considering the limit of constant m, in which case the term reduces to 
(/ pi'+m ) 2 } which is indeed the value of the depicted diagram. This procedure 
can be continued straightforwardly to obtain higher order diagrams. Let us 
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illustrate this for the second order term. In order not to make the formulas 
too lengthy we will denote the propagator by: 

G(x,y) =< x\G 1 \y > 
1 + mG 

The first and second derivative of the energy are now given by: 

f) P AT r 

— -r = -G(x, x) + 2N(N + 2)g J dyQ{y, y)Q(y, x)Q(x, y) + 0{g 2 ) (16) 

= d ^w^ + mN + 2)9 1 dyG{y -> xfG ^ y)2+ 

4N(N + 2)g J dyQiy, x)Q(x, y)Q{y, y)Q(x, x) + 0(g 2 ) (17) 

Again using this in (Q) we find after integration: 

AT r 
E= — Tr \og(l + mG) + N(N + 2) g J dxG(x,x) 2 + 

Ag 2 N(N + 2) J J dxd y g(y,x) 2 g(x,y) 2 + 

4g 2 N(N + 2) 2 J J dxdyQ(y, x)G(x, y)Q(y, y)Q(x, x) + 0(g 3 ) (18) 

Again by drawing the propagators or by considering the limit of constant m 
we find that this corresponds to the following diagramatic expression: 



4Af(iV + 2) 2 ff 2 CXX) + 4N(N + 2) g 2 




which are indeed the correct diagrams with the correct symmetry factors. 
It is clear that this procedure can be continued to any order one desires. Of 
course the expressions will become rather elaborate as one advances towards 
higher orders. For this reason we will present in the next section a com- 
pletely diagrammatic approach to the same problem, in which one does not 
work with mathematical expressions, but with Feynman diagrams that are 
combined together according to a set of simple rules. 
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3 Diagrammatic approach 



The most important formula in developping this approach is formula (]12|) for 
which we have already given a graphical expression. We will mainly repeat 
the steps followed in the preceding section, but this time from a diagrammatic 
point of view. The lowest order approximation of the energy is now given 
by: 



and the first functional derivative with respect to m(x) is given by (see for- 
mula (H)): 



dE N 



dm(x) 2 




+ 0(g) 



From now on we use formula (12). In the above diagram we have one prop- 
agator. Since the derivative doubles each propagator in turn (and changes the 
sign) we find the following diagrammatic expression for the second functional 
derivative: 



d 2 E _ N 



dm(x) 2 2 




+ 0(g) 



Taking the integral with respect to x in the diagrammatic approach cor- 
responds to identifying the different points that are labelled with " x" and 
dropping the label. By doing this we create a vertex. Doing so and using 
formula (ffl) we find: 

f = N(N + 2) 



recovering the result found in the preceding section. We can now proceed to 
higher orders by taking the functional derivative of the first order approxi- 
mation of the energy. This involves taking the derivative of the the two loop 
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diagram shown above. In this case we have two propagators, so the derivative 
will be the sum of diagrams, in which every line is doubled in turn (hereby 
changing the sign of the diagram). Diagrammatically, this gives: 

)=-*OQ- CO 



dm(x) 



x = — z V A 7 x 



The second order derivative is taken in the same way. We now have three 
propagators in each diagram, so doubling each of these in turn we arrive at 
the expression: 



a 2 



dm(x) 2 



) = 2 xiljx +4 



After integration with respect to x this gives: 
r dx—^ ( 

J dm(x) 2 \ 



X 



Using the bootstrap formula (f|) once more we find the diagrammatic equiv- 
alent of the expression found earlier in this paper: 

f = N(N + 2) OO + 



8N(N + 2fg QQQ+8N(N + 2)g +0{f) 

It is clear how this can be continued to higher orders. Both taking the 
functional derivative and integrating with respect to x are translated into 
simple diagrammatic manipulations, which allow a much simpler treatment 
than the one presented in the preceding section. Carrying on straightfor- 
wardly, one can easily generate Feynman diagrams to, in principle, an arbi- 
trary number of loops. We give here all diagrams up to five loops with the 
correct symmetry factor for general N: 
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+ N(N + 2)g 00 + 4iV ( iV + 2 )V 



+ 



4N(N + 2)g 2 (5) +16Ar(Ar + 2)VOOOO+647V(Ar + 2)vQDO 



+ 



f N(N + 2)(N + 8)g 3 + f N(N + 2fg 3 Q( 



64N(N + 2) 4 g 4 + 128 N(N + 4) V 




+ 



+ 



256iV(AT + 2)yOO(jD 

+ 256iV(A r + 2fg A 




38AN(N + 2) V 



+ 



+ 256iV(Ar + 2) 2 (A r + 8)# 4 QOO + 



32N(N + 2)VOOO + 32N(N + 2)(N 2 + 6N + 20)g 4 r ~, 

O 




64iV(iV + 2)(10iV + 44)£ 4 




+ 256N{N + 2) V 



4 Correlation functions 

So far we have shown how one can generate the diagrams that contribute to 
the energy of the model. If one were to use these diagrams to calculate the 
perturbation series, one would inevitably be faced with the problem of renor- 
malisation. In order to be able to carry out the renormalisation program, 
it is necessary to know the divergences of the two- and four-point functions. 
Therefore the technique presented in the preceding section would be of little 
use if it would not allow the generation of these diagrams also. For this rea- 
son (and for the sake of completeness) it is usefull to generalise the method 
to the generation of any kind of diagram (vacuum, two-point, four-point,...). 
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As a matter of fact this can be easily accomplished in a straightforward way 
by coupling a source term to the 0-field. The model is then defined by: 

1 ^ 
Z = e~ E = J £>0exp-- J dxd^d^cj) + m(x)(p 2 (x) - g J (<p 2 ) 2 - J J i & > 

■ — ■ ■ i=l 

(19) 

This extra term does not break the validity of (ffl), therefore the recursive 
procedure can be used without any change. The only change we have to deal 
with is the free field value of the energy. Instead of (^) we now find that 
E(g = 0) is given by: 

N 1 N 1 

E(9 = 0) = -Tr log(-9 2 + m) - - g J^ ^— J, , (20) 

where the "scalar product" in the last term is taken in coordinate space. 
Since the diagrammatic approach is far simpler and neater we will use it 
here from the very beginning. The above formula can be written graphically 

as: 



E(g = 0) = %() - \J J + 0(g) 



The same rules for derivation and integration over a point in space apply 
here. Taking the first and second derivative with respect to m(x) we find: 



8E N 



8m{x) 2 V»y 1 2 

X X 

X 

d 2 E N 



\J . J + 0(g) 



dm{x) 2 2 



- J ■ ■ J + 0(g) 

ry rp rp 



Putting this in (|j) and integrating with respect to g we find the first order 
value of the energy: 

J J 

E = E(g = 0) + (N(N + 2) OO +(iV + 2)J O J+j/\ j )g + Q(g 2 
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It is quite clear how the recursion can be continued towards higher orders. 
Doing this one generates E as a functional of the source J to an arbitrary 
number of loops. All correlation functions can then be obtained from this 
by taking successive functional derivatives with respect to the source. As an 
example we give here all one particle irreducible diagrams that contribute to 
the two point function up to four loops: 




Q_ + 



+ (N + 2)g* _0. 
ISi +3(A + 2)V- 



(JV+2)3 3 



2 




+ (A + 2)V 







+ (A + 2)V 





+ 



(JV+2) 4 „ 4 



2 



opp 



(N + 2) V 



+ 3(iV + 2)V 



+ 2(iV + 2)V 
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3(JV + 2)V ^-^ +3(JV + 2)V O +3{N + 2)Y 



o 




4(iV + 2) 2 (iV + 8)^ 4 ULi + (N + 2) 2 (N + 8)g i 
(iV + 2) 2 (iV + % 4 CjQ. + O 
(iV + 2)(iV 2 + 6iV + 20)# 4 C£Q + (iV + 2)(l(W + 44)c/ 4 
(AT + 2)(10iV + 44)(7 4 +6(iV + 2)V fjjN 



5 Gross Neveu model 

In this section we will present an analogous treatment of the Gross-Neveu 
model ||. This well known model is defined by the partition function: 

Z = e~ E = J ThpThjje-f^- 8 ^ (21) 

In order to obtain the perturbation series for this massless model we will 
consider the case where the field is allowed to have a mass m, which again 
depends on the coordinates. Instead of working with ([H]), we will investigate 
the theory with partition function: 

Z = e~ E = / D^e-J^^W^W) 2 (22) 



Once again we find a relation between the partial (functional) derivatives of 
the energy: 

— = [((—? - -^-) (23) 
dg J {K dm{x) ) dm(x) 2) 1 1 

The only missing ingredient is the value of the energy at zero coupling. Again 
without difficulty one finds that this is given by: 

E(g = 0) = -NTr\og(j8 + m) = -NTr log(l + mfl~ l ) + cte (24) 



12 



In this case "Tr" means the trace taken in coordinate space as well as the 
trace of the 7-matrices that occur in the expression (one could write it as: 
Tr = Tr^> ® Tr <yCOOrd \ where Tr^ denotes the trace of 7-matrices and 
rp r (coord) - g trace in coordinate space). This is an important difference 
with the first model considered in this paper. We can also derive the equiv- 
alent formulae of (|9|) and (|12"D, namely: 

OE 1 

iVTr (7) < x\G- -\x > (25) 



m{x) 1 + mG 

and: 

c) X XX 

< y\G- —\z >= - < y\G- — \x >< x\G- —\z > (26) 



dm(x) 1 + mG 1 + mG 1 + mG 



The last formula is in fact the same as (|T2]) , except for the fact that G is now 
JBT\ 

Instead of working with the rather cumbersome mathematical expressions, 
we will immediately turn to the diagrammatic approach proposed in the 
preceding section (with a few minor changes to make it better adapted to 
this model). The lowest order approximation of the energy could be written 
as: 



E = -N{) +0(g) 



It is however more convenient to attribute a minus sign and a factor of iV 
to the closed loop itself, so that the energy in diagrammatic notation would 
simply be: 

E = O 



We now proceed to generate the diagrams with this "initial value" and 
remember to assign to a certain diagram a power of iV equal to the number 
of loops in the diagram and a minus sign if this number is odd. Moreover, 
around every loop one would have to take a trace of the 7-matrices in the 
loop. Let us make this clear in the lowest orders. The first and second 
functional derivative are given by: 
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dm(x) 



X 



dm(x) 2 



d 2 E 



X 



o 



x +0(g) 



We still have the property that taking the functional derivative doubles a 
propagator in the diagram and changes the sign. The main difference is that 
the propagator is now the propagator of a fermion field. Using the bootstrap 
formula, we find for the energy of the Gross-Neveu model in first order: 



In the last formula we have slightly changed the convention used in the 
preceding section. Instead of merging two points with the same label, we now 
connect them with a dotted line, reminiscent of the auxiliary field notation. 
We do this in order to be able to keep track of the number of (index-)loops 
in each diagram, which is important for the "weight" of factors of N that 
the diagram will receive, as well as for the sign of the diagram. If one thinks 
of the dotted line as a delta-function it is easy to see that this is indeed 
the same as the convention in the preceding section. The next order in the 
coupling constant is found in much the same way: 



+(00+ ©)9 + eV) 




O 




x 



X 



and: 
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d 2 E 
dm(x)'- 



X 



+ 2g xO O x + 4g + 



x 



x 



■2n KZ) +0(g 2 ) 



which leads to the energy approximated in second order: 

% 2 O O +i6, 2 -vO:::::0 





The notation with the " auxiliary field" is clearly more convenient when it 
comes to taking traces of 7-matrices. The more conventional diagrams would 
obscure the way these traces have to be taken. When one wants to consider 
the massless case (as the Gross- Neveu model is originally defined) , a number 
of diagrams will automaticly yield zero. Indeed, in this case the propagator 
is given by Jr and therefore every propagator contributes one 7-matrix. This 
means that loops with an odd number of propagators automatically vanish, 
(this includes for instance all tadpole diagrams). Note that one cannot dis- 
card these diagrams before one has arrived at the end of the calculation. 
To obtain higher order diagrams correctly, one has to take into account all 
diagrams of lower order, also those that would vanish in the limit when m is 
taken to zero. Correlation functions for the GN-model can likewise be gen- 
erated by adding a source term. In the fermionic case however, one would 
have to introduce two source terms: 



Z = e 



-E 



= J VipVi/je~ S 



— j ip{^+m(x))il)— g(il)ip) 2 — r]il>— 1/17} 



(27) 



The lowest order approximation to the energy would then read: 
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E = Q -V V 



and one could use (^3|) to generate the energy as a functional of rj and rj. 
Correlation functions are easily found by taking successive derivations with 
respect to the source terms. 



6 Thirring model 

The third model which can be treated in a way similar to the vector model 
and the Gross-Neveu model is the Thirring model @. This model also is 
a theory of interacting fermions, but this time with an interaction term 
—gJ2^I{^l^) 2 - In order to be able to carry out the same techniques as for 
the previous model, it is not sufficient to consider an x-dependent mass. In 
this case we must take a mass that also depends on the Dirac indices. We 
will define this "matrix-mass" as: fn = r y fJ- m fJ- (x). The generalised partition 
function that we consider is now given by: 

Z = e~ E = J V~i)Vi)e- f^+^-g&i^) 2 ( 2 8) 

The formulas for the functional partial derivatives are given by: 

dE 



dm^x) 



and: 

sSr - < & ^ (x)f > + ■ (29) 

The bootstrap formula is easily generalised to: 

dE „ r , dE x2 d 2 E 

We also need to know the value of the energy when the interaction is switched 
of. Once more, this poses no serious problems and we find: 

E = -NTr log(/9+ jfn) = -NTr log(l + /rift' 1 ) + cte (31) 
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Just as in the Gross-Neveu model the trace is to be understood as a direct 
product of the trace in coordinate space and the trace of the 7-matrices in 
the expression. The formula for the functional derivative of an propagator 
is slightly different from the forms encountered above. We find without 
problems: 

dE \n 1 1 \n 1 1 \n 1 1 

d^Jx) < ylG T^G lz >= " < ylG T^G lx > ^ < xlG T^G lz > 

(32) 

This means that, apart from a doubling of the propagator and a changing 
of sign, the derivative inserts a 7-matrix at the point x. This poses no 
serious difficulties however in the diagrammatic expansion. As a matter of 
fact, one obtains the same diagrams as in the Gross- Neveu model (again 
using the notation with dotted lines to connect identical points). The only 
difference is the value that has to be assigned to each of these diagrams. 
The trace of 7-matrices around a loop now contains not only the matrices 
from the propagators in the loop, but also those inserted at the "vertices" of 
the fermion field and the " auxiliary field" . Because of this, the number of 7- 
matrices in a loop is always even. Therefore, unlike the massless Gross-Neveu 
model, all diagrams will contribute. The problem of generating diagrams that 
contribute to the various correlation functions can be treated by the addition 
of two extra source terms as in the case of the GN model. 



7 Non-local interactions - QED 

Another interesting generalisation of the method presented in this paper is 
the case of non-local interactions. Instead of an interaction term / dx(<f) 2 (x)) 2 
(or the analogon for the other models considered here) one could work with 
an interaction term: 

J dxdy<p\x)f\x,y)<p 2 (y) (33) 

where f(x,y) is some function of x and y. The diagrams that contribute 
to the energy of this model are not very much different from the diagrams 
of the models we have considered in the preceding sections. As a matter 
of fact, when we use the auxiliary field notation they are exactly the same! 
The only difference is that up to now the propagator of the auxiliary field 
was merely regarded as a notation for the ^-function and that in the end one 
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would have to identify the points that are connected by such a propagator. 
In the case of non- local interactions the propagator is a non-local function 
and the dotted lines now denote a real field propagating. This becomes clear 
when one considers the bootstrap equation for non-local interactions. We 
easily find: 

dE „/■,,,, w/ dE .. dE . d 2 E 

-=AJ dx dy /(s, y )((— )(_) - dm{x)dm{y) ) (34) 

If one uses the diagrammatic approach, one easily finds that the recursive 
method now comes down to marking two different points (x and y) in the 
diagrams and connecting them with a line which, in this case, represents a 
real field propagating between the two points. The functions f(x,y) can be 
seen as the propagator of the field in question. 

The most important example of a theory that has this kind of interaction 
is of course QED. QED is defined by the partition function: 

Z = e~ E = J VAViftfje" I \d^.dA»+w+m +91 »A^ (35) 

Since the action is quadratic in the photon field we can easily integrate this 
out and we end up with: 

g _ e -E _ J J)^p e - J ip(@+m)ip+f dx dyip-y l iip(x)A> J -'' (x-y)4>j„il>(x) j-gg^ 

where A is nothing but the photon propagator. This action is just a Thirring 
model with non-local interaction. Therefore the technique we have presented 
also allows one to generate Feynman diagrams that contribute to the energy 
and the various correlation functions of QED. 



8 Conclusion and outlook 

We have presented in this paper a method that allows to generate Feynman 
diagrams for theories that involve fields with N components, both bosonic 
and fermionic. By applying a simple algoritm one obtains diagrams of a 
certain order from the knowledge of the diagrams of lower order. A simple, 
purely diagrammatic, method for this can be constructed. In this paper we 
have restricted ourselves to quartic interactions, but one clearly sees that 
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interactions of higher order can also be treated in a similar way. In this case 
the equations @, and (|30"D will contain higher order derivatives with 
respect to m(x). For a interaction of order 2n one would have derivatives 
up to n th order. One is however always restricted to even order interaction 
terms (which was to be expected, since these are the only ones that make 
sense for general vector models with O(iV)-symmetry). 

The techniques presented in this paper apply only to vector models so far. 
It would of course be interesting to have a similar technique for the matrix 
models [|K|]. Other interesting extensions are the generation of diagrams 
in theories with odd interaction terms and theories with several interaction 
terms. The most important example of such a theory is of course QCD. 
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